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ABSTRACT 

Aims. Multiwavelength observations of the young massive cluster Westerlund 1 have revealed evidence for a large num- 
ber of OB supergiant and Wolf-Rayet binaries. However, in most cases these findings are based on the detection of 
secondary binary characteristics, such as hard X-ray emission and/or non-thermal radio spectra and hence provide 
little information on binary properties such as mass ratio and orbital period. To overcome this shortcoming we have 
initiated a long temporal baseline, multi-epoch radial velocity survey that will provide the first direct constraints on 
these parameters. 

Methods. VLT/FLAMES+GIRAFFE observations of Wdl were made on seven epochs from late- June to early- 
September 2008, covering ~35 confirmed members of Wdl and ~70 photometrically-selected candidate members. Each 
target was observed on a minimum of three epochs, with brighter cluster members observed on five (or, in a few cases, 
seven) occasions. Individual spectra cover the 8484-9001A range, and strong Paschen-series absorption lines are used 
to measure radial velocity changes in order to identify candidate binary systems for follow-up study. 
Results. This study presents first-epoch results from twenty of the most luminous supergiant stars in Wdl. Four new OB 
supergiant members of Wdl are identified, while statistically significant radial velocity changes are detected in ~60% 
of the targets. W43a is identified as a short-period binary, while W234 and the newly- identified cluster member W3003 
are probable binaries and W2a is a strong binary candidate. The cool hypergiants W243 and W265 display photospheric 
pulsations, while a number of early-mid B supergiants display significant radial velocity changes of ~15-25kms~ 1 that 
we cannot distinguish between orbital or photospheric motion in our initial short-baseline survey. When combined with 
existing observations, we find 30% of our sample to be binary (6/20) while additional candidate binaries support a 
binary fraction amongst Wdl supergiants in excess of ~40%, a figure that is likely to increase as further data become 
available. 

Key words, stars: evolution - supergiants - stars: binaries: general - techniques: radial velocities 
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| 1. Introduction radiation pressure can be overwhelmed, sufficient material 

g3 ■ must be pre sent to accrete from in order to yield the stars 

. . i Despite direct observational confirmatio n of stars with dy- themselves (|Davies et al Il2006h 
namical masses in ex cess of 80 M Q (e.g. iRauw et al. 1120051 : 

iSchnurret al.lfS) . the production of massive stars is Nevertheless, several differing scenarios have been ad- 
poorly understood, largely as a result of their intrinsic rar- vanced to surmount these difficulties. Disc mediated ac- 
ity, apparent rapidity of formation and extreme extinction cretion has been suggested to overcome radiation pressure 
due to veiling by their natal envelopes. Moreover, their pro- flYorke fc Sonnhalter | |2002J) but sufficient material must 
duction is problematic on theoretical grounds, primarily stl11 be supplied to the protostar. One mechanism pro- 
because accretion rates must be both extreme in order to P osed to accomplish this is competitive accretion ontohigh 
build the star and yet must overcome the resultant radia- mass Protostars in cluster cores (|Bonnell fc Bate | |2QQfl). 



tion pressure. Hence it is not clear if the process is simply a An alternative, which also circumvents the constraints un- 
sealed up version of low mass forma tion or whether it pro- P osed b y radiation pressure, is that massive stars form 
ceeds by a different mechanism (e.g. iBonnell fc BatellS as the result oOtellar mergers (e.g. | Bonnell et al. | I199J, 
and refs. therein). An additional problem is that even if Pavies et al. ||2006fl. though both processes may be rele- 
vant (|Bonnell fc Batell2005h . 

~S^d offprint requests to: B.W. Ritchie, e-mail: While direct observations of m site massive star forma- 

b.ritchie@open.ac.uk tion are technically highly challenging, significant observa- 

* This work is based on observations collected at the European tional constraints on this process exist. For instance, the 

Southern Observatory, Paranal Observatory under programme finding that mas sive s tars appear to form in clusters (e.g. 

ID ESO 81.D-0324A. . . E. IZinnecker et al. lll993tlClarke et al. 1 12000) appears difficult 
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to reconcile with turbulent fragmentation of the molecu- 
lar cloud leading to the production of protostellar cores, 
since the resultant clump masses are s ignificantly lower 
than required to form massive O stars (jClark fc Bonnelfl 
12004) . A second constraint is provided by the observed 
binary fraction and properties of massive stars; competi- 
tive accretion predicts their formation via accretion onto 
wide lo w mass systems leading to massive close binaries 
(-1 AU: iBonnell fc Bate l[2005fl . while massive star forma- 
tion via stellar mergers leads to a reduction in the number 
of primordial binaries, and hence requires a very high ini- 
tial fraction. As such, recent observational results demon- 
strating a large binary fraction amongst massive stars have 
the pote ntial to directly test current theories of star forma- 
tion (e.g. iGarcfa fc Mermilliod1 | 200l |. [Kiminki et al. Il2007l 
ISana et al. 1120081 iBosch et al. I bood IClark et al. Il2009ah . 
The properties of these populations, such as the distribu- 
tion of eccentricities, mass ratio and orbital separations 
place stringent constraints on the production rates and 
channels of both high and low mass X-ray binaries - and ul- 
timately double degenerate systems - as well as both Type 
Ib/c and blue Type II supernovae and po tentially Gamma 
Ray Bursters (jKobulnickv fc Fryer 1120071 and refs. therein). 

The starburst cluster Westerlund 1 (hereafter Wdl; 
IWesterlund I Il96lh contains a rich population of mas- 
sive, evolved stars including Wolf-Rayets, OB super- 
giants, yellow hy pergiants, a lumino us blue variable, and 
red supergiants (jClark et al. 1 l2005h . This unique popu- 
lation has been the subject of extensive study in re- 



Table 1. Summary of observations of the four sets of tar- 
gets. 



cent years at radio (iDoughertv fc Clark I 120081), infra- 



red (ICrowther et al. 1 120061). visual (IClark fc Neeueruela I 
120021 IClark et al 1 120051: iNeeueruela et al. I I2009T) a nd X- 
ray wavelengths (jSkinner et al. Il2006t IClark et al. Il2008t ). 
These studies have shown that 17 of the 24 i denti- 
fied Wolf-Rayet s tars are binary (ICrowther et al. [20061 : 



IClark et aDl2008D, while p hotometric (jBonanos 1120071 ) and 
X-ray (jClark et al. II2008T ) studies also imply a significant 



binary fraction amongst the transitional OB supergiants 
that are the evolutionary precursors of the Wolf-Rayet pop- 
ulation. However, the existing studies are largely based on 
secondary binary diagnostics such as the presence of non- 
thermal radio emission, hard and/or high luminosity X-ray 
emission and a significant near-mid IR excess, all of which 
are thought to be observational signatures of Colliding 
Wind Binaries (CWBs). Moreover, the single epoch obser- 
vations used to identify CWBs provide little information on 
the nature of the binary, save that the secondary is expected 
to be massive enough to support a powerful stellar wincQ. 
Binaries have been directly detected in only a few cas es, e.g. 
the eclipsing binaries W13 and W36 dBonanos 112007) or the 
doub le- lined spectroscopic binary W10 (jNegueruela et al. I 
12009ft . 

Consequently we have initiated a long temporal base- 
line, multi-epoch VLT/FLAMES spectroscopic radial ve- 
locity (RV) survey designed to identify and characterise bi- 
naries within Wdl. The sensitivity offered is high enough 
to identify both short- and long-period binaries and is sen- 
sitive to low mass companions within short period sys- 
tems. This is the first of three planned papers, describ- 
ing target selection, data reduction and analysis for our 



Date 


MJD" At 6 


Targets Integration time 


20/06/2008 
9Q /Dfi /9flflS 


54637.190 
54637.230 
e;zUSzUS rat; 

54646.224 — 


Faint 1 3 x 895s 
Faint 2 3 x 895s 

Rrio-ht 9 x fiflfk 

UIlgllL A A UUUS 

Faint 3 3 x 895s 


18/07/2008 


54665.036 18.85 
54665.065 18.84 
54665.104 27.87 
54665.144 27.95 


Bright 2 x 600s 
Faint 3 3 x 895s 
Faint 2 3 x 895s 
Faint 1 3 x 895s 


24/07/2008 


54671.134 6.10 


Bright 2 x 600s 


14/08/2008 
17/08/2008 


54692.042 20.91 
54692.073 27.01 
54695.094 29.99 
54695.133 29.99 


Bright 2 x 600s 
Faint 3 3 x 895s 
Faint 2 3 x 895s 
Faint 1 3 x 895s 


04/09/2008 


54713.011 20.97 


Bright 2 x 500s 



a The modified Julian day (MJD) is given at the midpoint of 
the sequence of integrations. 

b The time since the target field was last observed. 



survey, and presenting the first results from a subset of 
the VLT/FLAMES dataset containing a sample of lumi- 
nous, evolved supergiants observed on five epochs between 
29/06/2008 and 04/09/2008. Paper II will present analysis 
of our full datasetQ consisting of ~105 spectroscopically- 
and photometrically-selected targets observed on a mini- 
mum of six epochs during 2008 and 2009. A final paper 
will present follow-up observations a nd numerical model- 
ing o f the observed RV variability (c/. iKobulnickv fc Fryer I 
l2007t ) aimed at placing limits on the close, high-mass bi- 
nary fraction, the orbital and mass-ratio distributions and 
the permitted formation channels for low- and high-mass 
X-ray binaries within Wdl. 



2. Observations & data reduction 

Observations were made using the F ibre Large Array Mult i 



Element Spectrograph (FLAMES; iPasauim etld~f 12002). 



located on VLT UT2 Kueyen at Cerro Paranal. FLAMES 
provides multi-object spectroscopy using the GIRAFFE 
medium-high resolution and/or UVES high resolution spec- 
trographs, with the red arm of UVES fed by up to eight 
fibres and the MEDUSA mode of GIRAFFE allowing si- 
multaneous observation of up to 132 separate object^. 
Our observations use GIRAFFE in MEDUSA mode, but 
the high degree of crowding towards the core of Wdl 
means that we cannot simultaneously allocate fibres to all 
spectroscopically-confirmed cluster members due to limita- 
tions imposed by the physical size of the MEDUSA fibres 
and the need to avoid collisions. Instead we use FPOSS (the 
FLAMES Fibre Positioner Observation Support Software) 
to optimally allocate fibres to targets drawn from four lists 
of candidates: 



Wind-photosphere interaction may also result in X-ray emis- 
sion in systems where the secondary lacks a strong wind (e .g. 
the 09V + B1-1.5V binary CPD -41°7742: ISana et al. II2005T ). 



2 Additional data during 2009 is anticipated under ESO pro- 
gram ID 383.D-0633. 

3 Two GIRAFFE Integral Field Unit modes are also available. 
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Fig. 1. VLT/FORS1 R-band finding chart (~0.6" seeing) for the t argets listed in Ta ble H The field of view is ~4' square, 
corresponding to ~6pc at an assumed cluster distance of ~5kpc (|Clark et al. II2005T) . 



— Know n cluster members catalogued by IWesterlund I 
<|l987l ) . The brightest OB and transitional supergiant 
members of Wdl are assigned the highest priority for 
fibre allocation, while other catalogued OB I /ll stars 
are assigned moderate priority. Yellow hypergiants and 
the Wolf-Rayet population are assigned low priority; the 
YHGs are presumed too large to allow a close compan- 
ion, while the binary fractio n of the Wolf-Rayet pq pu- 
lation has been examined bv lCrowther et al. I (|2006f ). 

— Candidate cluster members with SuSI2 photometry 
(jClark et al. I |2005[ ). These targets are located in a 
5.5x5.5-arcminute field centred just north-east of the 
core of Wdl. These are again assigned moderate prior- 
ity, with increasing likelihood of selection away from the 
highest-priority OB supergiants in the cluster core. 

— Bright st ars away from the cluster core with WFI pho- 
tometry (|Baade et al. II1999T) suggestive of cluster mem- 
bership. These are assigned low priority, but still have 
a high probability of selection during fibre configura- 



tion due to the lack of alternative targets away from 
the cluster core. 
— Objects around the cluster core that are likely to be 
members based on their general colour, but have no ex- 
isting photometry or spectroscopy. These targets are as- 
signed a low priority for fibre configuration and serve as 
backup targets for fibres that cannot be allocated in any 
other way. 

Four lists of targets were generated using FPOSS. A bright 
list, observed on five epochs in 2008, contains 17 of the 
brightest confirmed cluster members and five candidate 
members located towards the core of Wdl for which no pho- 
tometry or spectroscopy exists. In addition, three faint lists, 
observed on three epochs, contain a further 17 known clus- 
ter members, 58 photometrically-selected candidates (28 
with SuSI2 photometry and 30 with WFI photometry) and 
16 candidates for which no existing observations are avail- 
able. A summary of all the observations is given in Table[TJ 
Each target list also includes ~10 fibres that were used 
to obtain sky spectra. Despite being assigned low priority, 
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Table 2. List of targets. 



ID 


kJLJCCLldl J- y IJvj 


RA ('I2000 , l 


V)pr ( T2000"! 


R" 


P 


Field 6 


Notes c 


W2a 


B21a 


i c Art cn 71 

lo 46 59.71 


-45 50 51.1 


14.23 


11.73 


Br. 


A 


W8b 


B1.51a 


~\ Ct /IT Cl a nr 

lb 47 04.95 


At- r n r\s* ry 

-45 50 26.7 






Br. 




W13 


B0.5Ia + +OB e 


16 47 06.45 


-45 50 26.0 


14.63 


12.06 


Br. 


X, E (9.20 day) 


W21 


B0.5Ia e 


16 47 01.10 


-45 51 13.6 


15.56 


12.74 


Br. 


A 


W23a 


B1.5Ia e 


16 47 02.57 


-45 51 08.7 


14.91 


12.07 


Br. 


A 


W30a 


O9-B0.5Ia bin / 


16 47 04.11 


-45 50 39.0 


15.80 


13.20 


Br. 


X 


W43a 


B01a e 


16 47 03.54 


-45 50 57.3 


15.22 


12.26 


Br. 


A 


W55 


BOIab 9 


16 46 58.40 


at r 1 1 1 n 

-45 51 31.2 


15.25 


12.67 


Br./l 


* 

A 


W60 


B01ab e 


16 47 04.13 


-45 51 52.1 


15.96 


13.28 


Br. 


A 

A 


W71 


B2.5Ia e 


-\ f Al~? no A A 

16 47 08.44 


A f r n < n o 

-45 50 49.3 


14.06 


11.16 


Br./l 


A 

A 


W78 


Blla e 


16 47 01.54 


A V A C\ V 7 O 

-45 49 57.8 


14.54 


12.04 


Br. 


A 

A 


W84 


(jy.oib 


Ad cn no 

16 46 59.03 


-45 50 28.2 


15.60 


13.63 


Br. 


A: 


W700 A 

W2J4 


uy.oib 


lb 4r 01.44 


-45 52 ,55.0 






Br./l 




W239 (F) 


WC9d 9 


16 47 05 21 


-45 52 25 


15.39 


12.90 


Br. 


X, A 


vv z,*±o 


T RV ( A ^Ta + , \' 1 
ijjD v i /\oia ( 


1 « 47 07 55 

ID *± i Ul .OO 


^9 98 ^1 






Rr 


A 


WzOO 


r ola 


-\a An c\« oa 
ID 4; uo.zD 


-40 4y zo. / 


lo.DZ 


in ^ 
1U.54 


D_ 

rsr. 


A 

/V 


Wo / o 


oUiaD 


1 a AG C7 71 
ID 40 / . / 1 


a k on n 
-40 Oo zU.U 


14. 14 




n_ 

r>r. 


A 


W3002 


BOIab 


16 46 54.24 


-45 51 54 7 






Br. 




W3003 


BOIb 


16 47 11.60 


-45 49 22.4 


16.21 


13.31 


Br. 


A 


W3004 


BOIab 


16 47 13.39 


-45 49 10.5 


15.96 


12.99 


Br. 


A 


W3005 


09.51b 


16 47 01.69 


-45 52 57.8 


15.39 


12.97 


Br. 


A 


Fl d 


M2II-III 


16 46 57.84 


-45 52 18.5 






Br. 


Field 


W6a 


O9-B0.5Ia bin/ 


16 47 03.04 


-45 50 23.6 


15.80 


13.16 


1 


X, P (2.20 day) 


W15 


09Ib e 


16 47 06.63 


-45 50 29.7 


16.38 


13.75 


2 


X, A 


W17 


09Iab e 


16 47 06.25 


-45 50 49.2 


16.19 


13.56 


1 


X 


W57a 


B4Ia e 


16 47 01.35 


-45 51 45.6 


13.83 


11.13 


1 


A 



"Photometr ic R and / -band photometric magnitudes are taken from Clar k et al. I (|2005T ) or iBonanosl (|2007T ) where available. 
Values from iBonanos (2007) are listed in italics. 

Targets are listed as Br. for members of the bright list, or 1. . . 3 for the faint target lists. The stars W6a, W15, W17 and W57a 

are not included in our main target list, but are discussed in t he text. 

c Stars are noted as Eclipsing, Periodic or Aperiodic variables (Bonanos 2007) and/or X-ray sources (|Clark et al. II2008T ). 

d The target Fl is a foreground star, and no t a member of Wdl. 

References: jNegueruela et al. I (|200g , 120091) ; ^ IClark et all (|2008l ); ^Clark et al.1 (I2005T ); ^Ritchie et all (|2009h . 



two Wolf-Rayets, W239(F/WR77n) and W39b(G/WR77j), 
were included by the fibre allocation process, with W239 
included in the bright targets list, while three targets are 
duplicated on the bright and faintl lists and are therefore 
observed seven times in total. With the exception of two of 
these duplicated objects (W71, B2.5Ia and W55, BOIab), 
the transitional supergiant W57a (B4Ia, faintl), and the 
Wolf-Rayet W39b/G (WN7o, faint2), the known cluster 
members in the faint li sts all have spectral clas sifications 
in the range 0911-BOIb (jNegueruela et al. Il2009t ), while vi- 
sual inspection of the photometrically-selected targets sug- 
gest they are predominantly reddened late-OI/II stars that 
are likely bona fide members of Wdl along with a number of 
field stars with colours similar to the reddened supergiants: 
full spectral classification will be given in Paper II. As many 
of our photome t rically -selected targets are uncatalogued 
by IWesterlund I (|l987t ) , we name targets that are subse- 
quently confirmed as members of Wdl as W2xxx if they 
are selected using photometry from lClark et al.l (|2005| ) and 
WSxxx otherwise, where in both cases the final three digits 
count upwards from 001. 

Due to the high de gree of reddening towards Wdl (A v ~ 
12: IClark et al. II2005T ) the expected signal-to-noise ratio at 
the blue end of the spectrum would be too poor to allow 
high-precision RV measurements to be obtained. Instead, 



we use the near-IR atmospheric window at ~ 8350 — 9000A, 
which covers a number of strong, well defined Paschen- 
series absorption lines in the spectra of O B supergiants 
(|Clark et al. Il2005t iNegueruela et al. II2009D . This re gion is 
free from telluric features, and the line formation region 
deep in the stellar photosphere implies the observed ab- 
sorption lines should be free from wind contamination^]. 
The GIRAFFE spectrograph was therefore used with setup 
HR21 covering the 8484-9001A range with R = A/AA ~ 
16200. Two 600s integrations were used for the bright 
targets, while three 895s integrations were used for the 
faint targets. The FLAMES data were bias-subtracted, flat- 
fielded to correct for pixel-to-pixel variation and fibre trans- 
mission differences, and wavelength-calibrated using ver- 
sion 2.5.3 of the FLAMES-GIRAFFE pipelinc0 with ver- 
sion 4.1.0 of the Common Pipeline Library (CPL), with 
individual spectra extracted from the final pipeline frames 



4 To test this hypothesis we have calculated non-LTE model 
atmosphere synthetic spectra for a range of parameters, with the 
synthetic spectra showing no evidence for wind contamination, 
even with mass loss rates an order of magnitude above expected 
values. 

http: / / www.eso.org/sci / data-processing/ software / pipelines / 
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Fig. 2. Spectral sequence of selected targets, from the M2II/III field star Fl to the 09.51b cluster member W84 with 
the rest wavelengths of the principal absorption lines marked. Residual sky lines, as discussed in Section [51 are visible in 
some spectra. 



using the IRAF^I task onedspec. The signal-to-noise ratio 
of our co-added spectra varies with the /-band luminosity 
of the target, ranging from S'/iV>200 for the very lumi- 
nous transitional hypergiants W243 and W265 to ~65-75 
amongst the least luminous 09.51b targets in our bright 
sample. A master sky spectrum was created from individ- 
ual fibre sky spectra, and this was subtracted from stellar 
spectra using the skysub task. Notable fibre-to-fibre varia- 
tions in the sky spectra are apparent across the FLAMES 
field, and as a result the removal of sky lines from stel- 
lar spectra is freq uently imperfect (see also discussion in 
lEvans et al. 1120051) . As can be seen from Figure [21 residual 
sky lines are not generally significant, but in a few cases 
unsubtracted sky lines overlapped absorption lines used for 
RV measurement and had to be removed by manually in- 
terpolating between levels on either side of the emission 
line. Finally, spectra were corrected for heliocentric veloc- 
ity using dopcor and normalized using the continuum task. 
A serendipitous check for zero-point errors in our data is 
provided by a strong, well-defined diffuse interstellar band 



6 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 



at -8620A (|M unari |[2000l ). As this feature is unassociated 
with Wdl it displays a constant profile and a RV that varies 
by less than ~lkms _1 in all spectra, implying that there 
are no systematic shifts in line centre or profile between 
epochs. 

In this work only the bright target list is examined; 
a list of targets is given in Table [5] and a finder chart is 
plotted in Figure [jQ. A representative spectral sequence, 
extending from the earliest objects in our bright sample 
(09.51b, e.g. W84 and W234) is plotted in Figure [H with 
the Paschen-series and He I lines labeled, along with the 
DIB at ~8620A, the C III A8500 line that blends with the 
Pa-16 line in stars of B0. 51 and earlier, and the neutral 
and singly-ionized metal lines that become prominent in 
the cooler stars. Radial velocities were determined from the 
strong Paschen-series and (when available) He I absorption 
lines, with the Pa-16 A8502 line excluded in stars earlier 
than BII due to blending with C III A8500 (this effect can 
be clearly seen in the bluewards shift in the Pa-16/C III 
blend in the 09.51b spectrum plotted in Figure [2]). In the 
case of the cool hypergiants W243 and W265, the Paschen- 



7 W373 lies outside the field of view of Figure \T\ and is located 
to the south of Wdl 16 arcseconds east of the bright field star 
HD151018. 
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series lines are frequently blended with adjacent Ca II and 
N I lines, and instead nine strong, unblended N I absorp- 
tion lines from high-excitation multiplets (Xlow ~10.2eV) 
we re used. Mea s ureme nts were made in a similar manner 
to iBosch et ai~l §009), by using the IRAF ngaussfit task 
within the stsdas package to fit Gaussian profiles to the 
absorption lines, with the measured stellar radial velocity 
at each epoch an error-weighted average of the individual 
absorption lines. Rest wavelengths for RV measurement are 
taken from the NIST Atomic Spectra databas^S Errors are 
<4kms _1 unless noted. 



3. Results 

The bright target list contains 16 stars that are known 
members of Wdl with prev i ously published spectral class i- 
fications (jClark et al. Il2005l.l2008t INegueruela et al. Il2008f) . 
The previously-unclassified target W373, located to the 
south of Wdl 16 arcseconds from the field star HD151018, 
is classified as BOIab based on the strength the Paschen- 
series and He I absorption lines, with a spectrum very 
similar to W60 (also BOIab; iNegueruela et aUl2009D . Four 
of the five candidate targets are confirmed as members of 
Wdl. W3002 and W3004 display very similar spectra to 
W373, and are again classified as BOIab. W3005 is slightly 
earlier, d isplaying a similar spec trum to W84 and W234 
(09.51b; INegueruela et al. I I2009D . while W3003 displays 
slightly broader Paschen-series lines than the BOIab stars 
and is classified as BOIb. Finally, although the fifth candi- 
date target (Fl) has a similar colour to the highly-reddened 
OB supergiants in Wdl, the lack of a well-defined DIB 
at 8620A and the strength of the TiO A8860 bandhead 
(|Ramsev 11198 If) show it to be a foreground early-M giant 
and not a cluster member. 

Table IA.1I lists the measured radial velocities for the 
target stars. In the case of W30a, the I-band spectrum 
shows only very broad, weak Paschen-series absorption lines 
that preclude accurate measurement; we return to this in 
Section [H Selected targets are examined further in the fol- 
lowing sections. Although the detection of binaries in Wdl 
is the long-term goal of our project, at this stage we are also 
interested in characterising possible photospheric sources of 
variations in RV. Both sources of variability are therefore 
discussed. 



3.1. The eclipsing binary W13 

iBonanosI (|2007h report W13 to be a contact eclipsing bi- 
nary with a 9.2-day orbital period, and it is also listed as a n 
X-ray source with L^-lO^ergs" 1 bv lClark et aTl (|2008ft , 
INegueruela et al. I (|2009h classify W13 as a B0.5Ia++OB 
binary, with one component being an emission-line object 
displaying strong Ha emission alongside C II AA6578, 6583 
and He I AA6678, 7065 lines in VLT/FLAMES LR6-mode 
spectra (jClark et al. I l2009bh . By chance, three of our ob- 
servations are separated by almost integer multiples of this 
orbital period: the second spectrum was obtained 2.05 or- 
bital periods after the first, while the fourth spectrum was 
taken after 4.98 orbital periods had elapsed. The other two 
are out of phase, taken after 2.71 and 7.26 orbital periods. 
We would therefore expect to see absorption lines at almost 
the same RV in three spectra with two discrepant, and this 
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Fig. 3. Pa-14 A8598 (all five epochs, left panel) and Pa- 
ll A8862 (three epochs at near-identical phase, right panel) 
in the spectrum of W13. The absorption feature at -8620A 
is a DIB, and rest wavelengths are marked. In this and 
subsequent figures the key gives the modified Julian date 
on which each spectrum was obtained. 



can be seen in Figure [3J An emission component is seen 
in all Paschen-series lines, weakening noticeably between 
the first and fourth spectra. Weak He I absorption features 
are also seen, with He I AA8583, 8777 (jvan Helden Ill972t ) 
most prominent. The radial velocities of the He I absorp- 
tion lines closely match the absorption components of the 
Paschen-series lines, and although no obvious He I emission 
is detected in our I-band spectra the expected strength is 
low. Radial velocity measurements and the virtual absence 
of an emission line component at phase ~0.7 (correspond- 
ing to the 24/07/2008, MJD=54671.1 spectrum) when the 
emission-line component is near (or at) eclipse show W13 
to be a spectroscopic binary, with the emission compo- 
nent originating in a B0.5I + /WNVL primary and the ab- 
sorption component in the secondary. The weakness of the 
Paschen-series absorption lines suggests a ^08-91 classifi- 
cation for the secondary, but this is inconsistent with the 
observed He I absorption lines which imply a later spec- 
tral type, and it is likely that there is significant infilling of 
the Paschen-series lines from the wind of the emission-line 
primary. W13 therefore appears to be an immediate pre- 
cursor to the binary-r ich Wolf-Rayet population in Wdl 
(|Crowther et al. 1120061 ). Radial velocity measurements and 
orbital parameters for W13 will appear in a forthcoming 
paper. 

3.2. The Wolf-Rayet W239 (F/WR77n) 

T he WC9d star W239 (F/W R7 7n in the nomenclatu res 
of IClark fc Negueruelall2002l and Ivan der Hucht 1 120061 re- 
spectively) has previously been highlighted as a likely bi- 
nary due to its strong near-IR excess and hard X-ray 
emission (ICrowther et aTl [2006. IClark et al.ll2008D . while 
iDoughertv fc Clark I ( 20081 ) report radio emission that is 
also consistent with a CWB. Our spectra are dominated 
by two strong C IIIAA8500,8664 emission lines, with weak 
C IVA8856 emission also apparent; the C III A8664 line and 
nearby DIB at ^8620A are plotted in Figured] Significant 
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Fig. 4. C III A8664 and the adjacent -8620A DIB in the 
Wolf-Rayet W239/F. A secondary component at A8653 re- 
sponsible for the extended wing bluewards of the main emis- 
sion peak is also marked. 




wavelength (A) 

Fig. 5. Pa-14 A8598 and adjacent ~8620A DIB in the spec- 
trum of W2a (top) and W43a (bottom). The Pa-14 rest 
wavelength marked. W234 and W3003, the two other can- 
didate binaries discussed in the text, are omitted for clarity 
and are intermediate between the two plotted spectra. 



radial velocity changes are apparent in all three emission 
lines, although the measured values of the two C III lines 
generally differ by ^20-40kms _1 , a result of excess emis- 
sion on the redwards side of the C III A8500 line shifting the 
profile fit relative to C III A8664 and also possibly the influ- 
ence of a weak secondary component at 8653A shifting the 
C III A8664 fit bluewards. Measurement is also complicated 
by the emission line profiles, which are sometimes skewed 
or weakly double-peaked. However, over the course of our 
observations RVs measured from the C III A8664 line span 
a range of -82kms _1 to +14kms~ 1 , with the C III A8500 
line giving a very similar overall range that is redshifted by 
~25kms _1 (errors estimated at <10kms _1 in both cases), 
with a large delta of ~65kms _1 between the 24/07/2008 
(MJD=54671.1) and 14/08/2008 (MJD=54692.0) spectra. 
Although WRs typically display low-level variability in 
emission lines due to clumping in the wind, the RV changes 
seen here are strongly suggestive of binarity. Dust-forming 
WC stars are generally con sidered to be binary systems 
with an OB companion (e.g. iTuthill et al.~lll999h . and the 
observed RVs are cons istent with a ~15M Q WR primary 
(jCrowther et al. Il2006| ) and more massive OB secondary in 
a short-period orbit viewed at moderate inclin ation (un- 
like B /WR77o, eclipses are not visible in W239: IBonanos I 
2007). The limited sampling of our five epochs precludes 
estimation of an orbital period, but our long-term dataset 
will allow us to constrain the orbital parameters of W239. 

3.3. The spectroscopic (candidate) binaries W2a, W43a, 
W234 and W3003 

Figure [5] plots the Pa-14 A8598 line and the adjacent DIB 
at 8620A for W2a and W43a, with both objects displaying 
notable RV variations in the photospheric absorption lines. 
W43a is unambiguously a short-period binary, with the 
measured radial velocities spanning a range of ^140kms _1 
and including a delta of ^131kms _1 in the space of six 
days between 18/07/2008 (MJD=54665.0) and 2 4/07/2008 

MJD=54671.1). No eclipses are reported by IBonanos! 

2007) , but the relatively high radial velocities suggest that 



the system must be near the sin i~0.7 limit where eclipses 
would become visible. The ~35M Q initial mass of the BOIa 
primary requires a massive companion to produce the ob- 
served RV changes, and although there is not direct spec- 
troscopic evidence for this secondary, changes in the He I 
line profiles are tentatively suggestive of an OB compan- 
ion. However, we note that changes in the Paschen-series 
absorption line strengths are also apparent in our spectra, 
and it seems likely that we are also observing photospheric 
variations superimposed on the changes in line centre due 
to orbital motion. This may explain changes in the He I pro- 
file, and follow-up observations at other wavelengths may 
be required to identify the companion. Nevertheless, W43a 
represents an encouraging prospect for accurate parameter 
estimation as further data become available. 

W234 is duplicated on the bright and faintl lists and 
therefore observed on seven occasions. The overall range 
of observed RVs is far smaller than W43a (~36kms _1 ), 
but relatively rapid RV changes are again present with 
spectra from 18/07/2008 (MJD=54665.0) and 24/07/2008 
(MJD=54671.1) showing a delta of ~35kms- 1 . This is 
large for bulk photospheric motions in the Paschen-series 
line forming region, and it is more likely that W234 is a 
similar short-period binary to W43a but viewed at a less 
favourable angle for RV measurement. W2a displays RVs 
ranging from -SG^kms" 1 to -18.6kms _1 over the course of 
our observations, including a change of ~22kms~ 1 within 
the six days separating the 18/07/2008 (MJD=54665.0) 
and 24/07/2008 (MJD=54671.1) spectra, and a change of 
^Bkms^ 1 between the final two spectra. The variability 
on short timescales is less pronounced than in W234, but 
the amplitude of the RV changes is still large relative to the 
majority of our sample. However, we cannot categorically 
rule out a photospheric origin for these variations and clas- 
sify W234 as a probable binary based on the magnitude of 
its short-term RV changes and W2a as a candidate binary 
pending further observation. 

Finally, the newly-confirmed cluster member W3003 
displays RVs spanning ^50kms _1 over the course of our 
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3.4. The pulsating hypergiants W243 and W265 
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Fig. 6. Spectral variability in W265. The left panel plots 
the region around the Pa-14 A8598 line, while the right 
panel plots the region around the N I triplet at A8680. . . 6. 
Rest wavelengths are indicated. 
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Fig. 7. Spectral variability in W265. The left panel plots 
Fe I A8824 (multiplet 60) and A8828 (multiplet 1269), with 
velocities relative to the Fe I A8824 rest wavelength. The 
right panel plots the Ca II A8498/Pa-16 A8503 blend, with 
velocities relative to the Ca II A8498 rest wavelength. 



observations. Unlike the other short-period (candidate) bi- 
naries identified here, W3003 shows no rapid deltas in 
RV but increases from -4.4kms _1 in the first spectrum 
(29/06/2008, MJD=54646.2) to -54.6kms- 1 in the third 
(24/07/2008, MJD=54671.1) before turning around and de- 
creasing to -10.6kms _1 in the final spectrum (04/09/2008, 
MJD=54713.0). Like W234, additional observations are re- 
quired to confirm the binary nature of the system, but the 
pattern of RV changes appears best explained by a limited 
sampling of a longer-period orbit viewed at moderate incli- 
nation. We therefore classify W3003 as a second probable 
binary. 



W243 (LBV ~A3Ia +: iRitchie et al. I I2009T) and W265 
(F5Ia+; IClark et al. I l2005h represent two of the popu- 
lation of cool hypergiants in Wdl. Turning first to the 
yellow hypergiant W265, the VLT/FLAMES spectra re- 
veal pronounced spectral variability in the neutral metal 
lines, and in particular in low-excitation multiplets of Fe I, 
Mg I and Si I. This can be seen in Figure [HJ which 
plots the regions around the Pa-14 A8598 line and the 
N I triplet at ~8683A on two epochs, separated by ^46 
days, that correspond approximately to apparent mini- 
mum and maximum T cS (29/06/2008, MJD=54646.2 and 
14/08/2008, MJD=54692.0 respectively). The first spec- 
trum displays strong Fe I lines from multiplets 60 (a 5 P- 
z 5 P°, Xiow ~2.2eV), 339 (b 3 P-z 3 P°, xiow ~2.8eV) and 401 
(b 3 G-z 3 G°, xiow ~3eV), but in the later spectrum the mul- 
tiplet 60 line (A8689) has weakened substantially while the 
multiplet 339 lines (AA8611, 8675) are undetectable and 
the multiplet 401 line (A8582) is weakly in emission. This 
pattern continues in the left panel of Figure Ofl which 
plots the adjacent Fe I A8824 (multiplet 60) and A8828 
(multiplet 1269, x 5 D°-e 3 D, xiow ~4.9eV) lines over all five 
VLT/FLAMES epochs. Once again, the multiplet 60 line is 
seen to vary significantly in both strength and profile, while 
the multiplet 1269 line strengthens in emission. In contrast, 
the N I, Ca II and Pa-16 lines visible in Figures [6] and [7] 
display much less variability. The Fe I lines (xiow ~2-3eV) 
and Mg I lines (xiow ~4eV) form in the upper photosphere, 
while the N I (xiow ~10eV) and Paschen-series form at 
lower levels. Notable differences between the two layers are 
apparent in Figure [5J which plots the measured radial ve- 
locities of two low-excitation and two high-excitation lines 
over the course of our observations. The outer photosphere, 
traced by low-excitation Fe I and Mg I lines, displays a 
pronounced infall followed by rapid expansion in the fi- 
nal epoch. The observed minima in the strengths of the 
Fe I absorption lines corresponds to maximum infall (and 
presumably maximum T c g), with these lines subsequently 
broadening due to the rapid expansion in the final epoch. In 
contrast the lower layers of the photosphere are more sta- 
ble, and although all species display similar radial velocities 
during the initial infall phase, the high-excitation lines show 
contraction ending earlier with a less pronounced expansion 
in the final epoch. 

The behaviour of the Fe I lines is remarkably simila r 
to the well-studied YHG p Cassiopeia (jLobel et al. I [20031. 
with the multiplet 60 line in Figure [7] developing a sim- 
ilar triangular profile wit h excess absorption in the blue 
wingTJ3 to that reported bv lLobel et al. I (|l998h for the low- 
excitation Fe I A5572.8 line. Our limited spectral coverage 
prevents us f rom examining othe r low-excitation absorption 
lines used bv lLobel et al. I (|1998h to interpret the spectrum 
of p Cas, but, given the similarities, it is likely the Fe I 
lines share a common origin in a (non-radially) pulsating 
photosphere. 



Note that for clarity only four epochs of data are shown in 
the right panel, with the omitted 17/08/2008 (MJD=54665.0) 
profile very similar to the plotted profile from 24/08/2008 
(MJD=54671.1). 

10 The excess blue-wing absorption forms in a co ol, optically- 
thick expanding wind at smaller optical depth ijLobel et al. I 
l200l l. 



B. W. Ritchie et al.: Massive binaries in Wdl 



9 




54670 54680 
Modifi9d Julian Day 



Fig. 8. Heliocentric radial velocity changes in low- and 
high-excitation absorption lines in W265. Plotted are 
Fe I (multiplet 60, xiow ~2.2eV), Mg I (multiplet 7, 
Xiow ~4.3eV), N I (multiplet 1, xiow ~lQ.2eV) and the 
Paschen-11 line (xiow=12.04eV). Errors are ~ ±2kms . 
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Fig. 9. Spectral variability in the LBV W243. The left 
panel plots N I A8629 line and adjacent ~8620A DIB from 
three epochs (with the other two omitted for clarity) while 
the right panel plots the N I A8712 line from all five epochs. 
Rest wavelengths are marked. 



Finally, we note the weak emission evident in the core 
of the strong Ca II A8498 line in Figure [3 a similar effect is 
seen in the core of the other Ca II AA8542,8662 multiplet-2 
lines, but in these cases infilling of the line centre redwards 
of the emission peak is seen, leading to an asymmetric core 
profile. The 4 2 P° (xiow ~1.7eV) lower level of the near-IR 
Ca II triplet is fed by the Ca II-K and -H lines in the near- 
UV, and core emission in the Fraunhoffer lines is an indica- 
tor of chromospheric activity in cool stars. Such cores are 
reporte d in a Orionis but are not observed in the near-IR 
triplet (jLobel fe Dupree Il2000f) . implying a higher chromo- 
spheric temperature in W265. Observations of p Cas show 
that core emissi on in the chromosp heric Ca II-H and K lines 
is not constant (jLobel et al. II2003T ). implying that a perma- 
nent chromosphere is not present, and this may also be the 
case for W265 ; no apparent chro mospheric X-ray emission 
was detected (|Oark et al. II2008T section 4.2.3), suggesting 
that the chromospheric activity implied by the Ca II lines 
may also be transitory. Further observation and modeling of 
W265 is required to determine the origin of these features. 
lLobel et al. I (|2003[ ) report line-splitting in low-excitation 
(Xiow <1.5eV) Fe I multiplets and other low-excitation met- 
als. However, our limited spectral coverage precludes exam- 
ination o f the m ultiplets with in which splitting is reported. 
iSargent I (|1961h lists Fe I multiplet 60 as containing dou- 
ble lines in p Cas, but no evidence of core doubling in this 
multiplet is seen in our spectra of W265. 

Significant pulsational variability is also observed in 
W243. The highly- variable Fe I lines observed in W265 are 
not present in the hotter LBV (A3Ia + ), which displays a 
complex Fe II spectrum along with absorption lines from 
other neutral and singly-ionized metals and an emission- 
line spectrum containing Balmer- and Paschen-series lines, 
He I and Ca II. This emission-line spectrum is a result of a 
hot companion star ionizing the wind of the LBV primary 
(|Ritchie et al. |[2~009f ). but this is not directly reflected in 
the cool-state absorption lines which display radial veloci- 
ties in a narrow range of -45kms _1 to -lSkms^ 1 across ten 
epochs of VLT/UVES and VLT/FLAMES data. These lines 



display RV variations on timescales of days, but - if they 
were due to orbital motion - we would expect a much wider 
range of radial velocities than observed unless, by chance, 
we are observing the system from almost directly above 
the orbital plane. The variations are most clearly seen in 
the Si II AA6349, 6371 doublet (obtained using VLT/UVES 
with cross-disperser #3, R~40,000) and in the near-IR N I 
lines; an example of the latter is shown in Figure In ad- 
dition, VLT/UVES observations of the N I lines reveal the 
formation of excess bluewards absorption similar to that 
seen in the Fe I multiplet 60 line in W265 at a time when 
Fe II emission lines develop pronounced P-cygni profiles in- 
dicating enhanced mass loss. W243 therefore appears to be 
undergoing very similar pulsational mass- loss to W265, al- 
beit at higher T e g; we exami ne the spectrum of W243 in 
detail in lRitchie et al. I (|2Q09h . 

3.5. The early-B supergiants 

Statistically-significant RV changes are also apparent in 
a number of the early-B supergiants, including W8b, 
W21, W23a, W71 and W3004. Of these, W71 (B2.5Ia; 
iNegueruela et aL~1l2009T ) is one of the latest of the continu- 
ous sequence of OB supergiants in Wdl, displaying strong 
Paschen-series and He I absorption lines along with weak 
N I lines indicative of its cooler state. W8b and W23a are 
slightly earlier, displaying no apparent N I lines. Finally, 
W21 and W3004 are the earliest of this group, with spec- 
tral types B0.5Ia and BOIab respectively. With the excep- 
tion o f W8b all are listed as aperiodic variables bv lBonanos I 
()2007f ). and none are sig nificant detections at X-ray wave- 
lengths (jClark et al. Il2008f) . 

In all cases, the observed RV variations lie in a narrow 
range of ~15-25kms _1 with no rapid variations between 
epochs. Although this is compatible with longer-period or- 
bital motion, the radial velocity changes in the Paschen- 
series lines are accompanied by changes in absorption line 
strength and profile, and it is more likely that the origin of 
the RV changes is photospheric. This is best illustrated by 
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Fig. 10. The Pa-14 A8598 line and adjacent -8620A DIB 
in the spectrum of W57a (top) and W71 (bottom). The rest 
wavelength of Pa-14 is marked. 
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Fig. 12. As Figure [TT} but showing variations in the 
Paschen-series absorption lines in the short-period binary 
W43a. 
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Fig. 11. Variations in the Paschen-series absorption lines 
in W78. The left panel plots the Pa-14 A8598 line and ad- 
jacent He I A8583 line, with rest wavelengths marked, while 
the right panel plots the Pa- 11 A8862 line with velocities 
relative to its rest wavelength. 



W71 and the cooler tr ansitional supergiant W57a (B4Ia; 
iNegueruela et aL~1l2009f ) which is included in the faintl list 
and was observed on three occasions; both are plotted in 
Figure [TDJ As well as variations in the Paschen-series lines, 
both objects also show pronounced changes in the strength 
of the N I A8680 triplet over the course of our observa- 
tions, a sensitive indicator of temperature that streng thens 
rapidly from B2Ia to B4Ia (jNegueruela et al. 1 12009) and 
is therefore strongly suggestive of photospheric changes. 
In addition, W78 (BHa) does not display statistically sig- 
nificant RV changes (<10kms _1 during the course of our 
observations) but nevertheless displays notable changes in 
the strength of the Paschen-series absorption lines. These 
changes, plotted in Figure ITTT are not apparent in the ad- 
jacent He I line which again suggests a photospheric origin. 
Radial velocity or line profile changes are therefore observed 
in all stars later than ^B0.5Ia in our sample, and it is possi- 



ble that the early- to mid-B supergiants in our sample all lie 
within a region of pulsational instability that extends to th e 
cool hypergiants (|Maeder fc Rufener II i 9 72t I Schaller 1l 1 9 900 . 
Unless the stars are short-period binaries at favourable in- 
clination (as with W43a, which displays clear signs of both 
orbital and pulsational variations; see Figure fT2|) we can- 
not conclusively detect binarity with our short baseline of 
observations, although the longer baseline survey will be 
sensitive to long-term trends in these objects. 



4. Discussion and Conclusions 

4.1. Binary detections 

In this work we present first-epoch results from a 
VLT/FLAMES radial velocity survey designed to identify 
high-mass binaries in the starburst cluster Westerlund 1. 
As expected, large RV change s are detected i n the short- 
period eclipsing binary W13 (|Bonanos I l2007h which dis- 
plays variable emission from a B0.5Ia + /WNVL emission 
line primary and He I and infilled Paschen-series absorp- 
tion in an early-B secondary. As such, it appears to repre- 
sent the immediate evolutionary precursor to the binary- 
rich Wolf-Rayet p opulation in Wd l, such as the eclipsing 
binary B/WR77o (|Bonanos Il200l . The WC9d star W239 
(F/WR77n) also shows significant RV changes compatible 
with a short-period WR+OB CWB viewed at moderate in- 
clination: the OB seconda ry is already strongly i mplied by 
the presence of hot dust (jCrowther et al. 1120061) , hard X- 
ray emission (also seen i n W13: | Clark et al. Il2008fl a nd non- 



thermal radio emission ([Dougherty fc Clark Il2008h . In ad- 
dition we identify W43a as a new short-period binary with 
a massive, so-far unresolved companion. W234 and W3003 
are probable binaries, and W2a is a strong candidate bi- 
nary. Both W2a and W234 appear to be short-period bina- 
ries, with W234, at spectral type 09.51b, earlier than the 
apparently pulsationally-variable early-B supergiants dis- 
cussed in Section [331 Finally, the measured RVs for W3003 
are consistent with limited sampling of a short- or longer- 
period binary. Future observations will resolve the nature 
of these systems. 
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An second important group of objects in the context 
of our RV survey are W6a, W15, W17 (included in the 
faint lists) and W30a, which have all been identified as 
strong binary candidates in previous studies of Wdl. W30a 
is the strongest X-ray so urce in Wdl after the magnetar 
CXOU J164710.2-455216 (jClark et al. II2008T) and a factor 
of ^10 more luminous than would be expe cted from a sin - 
gle staiPI. while W6a is a periodic variable (|Bonanos 1120071 ) 
that is also a hard X-ray source. Both display almost fea- 
tureless I-band spectra with very broad, shallow Paschen- 
seri es absorption lines typica l of spectroscopic binaries (see 
also iNegueruela et al. 1120091 ) and broad, variable Ha emis- 



sion (IClark et al.J|2008l ). W15 and W17 are both detected 



at radio ( Dougherty fe Clark ||20 08) wavelengths: the spec- 
tral index of W15 suggests composite stellar+non-thermal 
emissiorQ, while the index of W17 is unambiguously non- 
thermal. The radio observations therefore suggest that both 
are CWBs, altho ugh they are only weak detections at X- 
ray wavelengths ([Clark et al. Il2008h . None of these objects 
show significant RV variations in our data. In the case of 
W6 and W30a this is unsurprising, as the broad, weak I- 
band spectral features are insensitive to changes in RV due 
to the difficulty in determining the line centre, although 
at favourable inclination the large RV changes expected of 
a short-period, high-mass binary would still be apparent. 
In contrast, W15 and W17 display well-defined Paschen- 
series absorption lines but no significant variations in RV. 
However, in these cases it is likely that wider orbits are 
required in order for non-thermal radio emission to be de- 
tectable, and the relatively short baseline of our first obser- 
vations remains insensitive to longer-period, radio-strong 
binaries. 



4.2. Pulsations and the cluster velocity dispersion 

The cool hypergiants W243 and W265 are both undergoing 
photospheric pulsations, with W265 showing strong spec- 
troscopic similarities with the non-radia lly pulsating yel- 
low hypergiant p Cas (jLobel et al. 1 120031 ) . Chromospheric 
activity is also apparent from core emission in the near- 
IR Ca II triplet. We also find significant RV variations 
of ^15-25kms _1 in five early-B targets: while these could 
be orbital, in all targets later than spectral type ^Blla 
these RV changes are accompanied by changes in absorp- 
tion line profile suggestive of photospheric pulsations. The 
two-month baseline of our initial data offers prevents us 
from distinguishing a limited sampling of a long-period 
orbit from photospheric pulsations in the early-B super- 
giants, and further observations are required to break this 
degeneracy. However, regardless of their origin, these vari- 
ations imply that care must be taken in deriving a stellar 
velocity disp e rsion for Wdl from a single epoch of data. 
iBosch et al. I ()2009f ) find that the binary population leads 
to a significant overestimation of the velocity dispersion 
for NGC 2070 unless binaries are removed from the sam- 
ple, and the pulsational variations reported here may have 
a similar effect on infra- red determinations of the velocit y 
dispersion for Wdl fe.g. lMengel feTacconi-Garman | [2008). 
We note that separate VLT/FLAMES observations of the 
F2Ia + YHG W4 also reveal apparent pulsational variations 
in the Fe I lines (jClark et al. H2009bT ) and consider it likely 
that photospheric pulsations are a general characteristic of 
the transitional B-type supergiant and cool hypergiant pop- 
ulations of Wdl. 



Radial velocity surveys are also expected to be incom- 
plete, as a subset of objects will be viewed at unsuitable 
orbital inclinations while others do not display spectra that 
allow accurate RV measurement. There is also an element 
of chance in conclusively detecting binarity from a lim- 
ited number of RV measurements, as is apparent from the 
observations of W13 and W43a: only one of five observa- 
tions W43a unambiguously identifies it as binary, while by 
chance three observations of W13 fall at almost exactly 
the same orbital phase. However, while the short-per i od bi- 
nary W13 was detected photometrically (|Bonanos I |2007[ ) 
and wa s also as a strong binar y candidate f rom spectro- 
scopi c (jNegueruela et al.l l2009t ) and X-ray (jClark et al. I 
2008) surveys, neither W43a or the probable/candidate bi- 
naries W2a, W234 and W3003 have been identified in pre- 
vious work. Therefore, although the principle aim of our 
VLT/FLAMES survey is to build a dataset that is sensitive 
to long-period binaries, there is also great value to close- 
spaced RV measurements (e.g. the pair from 18/07/2008, 
MJD=54665.0 and 24/07/2008, MJD=54671.1, separated 
by six days) that effectively detect short-period binaries 
that are not apparent from photometry or single-epoch 
spectroscopy. 



11 The optica l spectrum of W30 a is also indicative of an inter- 
acting binary ([Clark et al. 1120081 ). 

12 Purely thermal emission from W15 would imply a very high 
mass loss rate inconsistent with the observed spectral type. 



4.3. Binary fraction 

Despite the expected limitations of our short-baseline ini- 
tial dataset, the results presented here support the high 
binary fraction implied from the earlier photometric and 
multiwavelength observations of Wdl. Considering just the 
sample presented here, RV measurements show W13, W239 
and W43a are binaries with W234, W3003 and W2a likely 
candidates. Radio and X-ray measurements strongly sug- 
gest a further four targets (W6, W15, W17 and W30a) 
are binaries but RV variations are undetected - and, in 
two cases, potentially undetectable - in our VLT/FLAMES 
data. Any signature of orbital motion in the two cool hy- 
pergiants is masked by photospheric pulsations, but the 
LBV W243 requires a hot companion to produce the ob- 
served He I and Lya-pumped Fe II and O I emission lines 
that are incompatib le with the cool-phase LBV primary 
(jRitchie et al.ll2009D . The YHG W265 is also plausibly 
binary, being associated with a compact, resolved radio 
source ([Dougherty fe Clark II2008T ) which could not be ion- 
ized by the YHG itself. While the source of ionizing pho- 
tons could simply be the OB supergiant population of Wdl, 
the YHG HR 8752 is also detected at radio wavelengths 
as a result of the wind being ionized by a B1V co mpan- 
ion ( Stickland fe Harmer I [l978t iPiters et al.lll988D . The 
hot companion also explains the observed [N II] emission 
in the spectrum of HR 8752; notably [N II] AA6548, 6583 
emission is visible in VLT/FLAME S LR6-mode spect ra of 
W265 (but not in the other YHGs: IClark et al. Il2009bl) and 
[N II] 
(fl987l) 



A6583 emission was also reported by IWesterlund I 
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Therefore even a conservative estimate that includes 
just the robust RV detection^] and the two objects that 
lack RV detections but have binarity strongly implied by 
other observation^ suggests a binary fraction of 30% 
(6/20) amongst our initial bright sample of supergiant stars 
in Wdl, increasing to >40% if other candidate binaries are 
included^. Both percentages are likely to increase as more 
data become available. If the more evolved WR/WNVL bi- 
naries and transitional hypergiants are excluded, then we 
find ~19% (3/16) of the OB supergiants in our sample are 
strong binary detections, rising to ~30% (5/16) if candi- 
dates are included. Robust estimates of the binary fraction 
amongst the OB supergiant population of Wdl will be- 
come available once the full FLAMES dataset is analysed 
in Paper II of this scries. However, these values are broadly 
consistent with the high binary fr action in the Wolf-Rayet 
population (|Crowther et al. |[2006D ; future observations and 
Monte-Carlo simulations planned for Paper III will reveal 
if the binary fraction of the two populations are equivalent. 
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Appendix A: Radial Velocity measurements 



B. W. Ritchie et al.: Massive binaries in Wdl 
Table A.l. Error-weighted radial velocity measurements for the targets listed in Table ^ a,h 
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ID 


Epochs 


Line(s) c 


Range d 
(kms" 1 ) 


54637 


54646 


Modified Julian Date 
54665 54671 54692 


54695 


54713 


Notes 


Wza 


r 



Fall-lo, Me 1 


Kin i/ii 
— ol.U 14.1 




A CI O 

— 4y. 6 


-56.2 


-33.9 


-54.1 




i o a 
— 18. D 


Cand. binary 


V V O D 





Do 11 1 c 

.rall-lo 


— Dl.D 4U.O 




a i a 
— Dl.D 


-49.9 


-40.3 


-47.4 




— OU. / 




W91 


c 



Pc»1 11/1 


— ( — ol.l 




— t .*± 


-35.0 


-43.1 


-55.1 




*}1 1 
— Ol.l 




VV ZiOoj 


o 


Pall-16 






—47 Q 


-55.6 


-56.6 


-42.9 




—62.4 




W43a 


5 


Pall-14 


— 1 37 Q — 1-8 1 




— 109.7 


-123.1 


+8.1 


-38.8 




— 137.9 




W55 e 


7 


Pall-14 




—43.1 


—49.9 


-55.4 


-45.4 


-43.8 


—53.9 


—55.5 




W60 e 


5 


Pall-14 






—40.5 


-48.3 


-45.8 


-48.2 




—37.7 




W71 


7 


Pall-16 He T 


_fil q S8 7 


—57.3 


—52.9 


—39.4 


—47.3 


—44.8 


—46.4 


—62.4 




W78 


5 


Pall-15 






—46.7 


-46.4 


-41.3 


-39.9 




—49.5 


Puis RSO 


W84 


5 


Pall-13 






-53.8 


-59.0 


-54.9 


-54.6 




-55.8 




W234 e 


7 


Pall-13 


-67.0 - -30.8 


-67.0 


-66.9 


-66.2 


-30.8 


-46.4 


-56.9 


-32.7 


Prob. binary 


W239 / 


5 


C III A8664 


-81.6 - +14.1 




-26.9 


-10.7 


+14.1 


-49.7 




-81.6 


Wolf-Rayet 


W243 9 


5 


N I (1 & 8) 


-28.8 - -20.4 




-22.7 


-20.4 


-27.7 


-28.3 




-28.8 


Puis. LBV 


W265 9 


5 


N I (1 & 8) 


N/A h 




-47.6 


-36.8 


-36.0 


-37.9 




-49.8 


Puis. YHG 


W373 e 


5 


Pall-14 






-50.8 


-50.9 


-45.8 


-53.8 




-52.9 




W3002 e 


5 


Pall-14 






-36.1 


-38.6 


-47.7 


-39.4 




-43.0 




W3003 e 


5 


Pall-14 


-54.6 - -4.4 




-4.4 


-43.7 


-54.6 


-22.0 




-10.6 


Prob. binary 


W3004 e 


5 


Pall-14 


-62.3 - -35.5 




-47.9 


-62.3 


-54.3 


-35.5 




-46.0 




W3005 e 


5 


Pall-13 






-43.9 


-46.9 


-52.4 


-51.6 




-49.9 





"Errors are not listed, but are <4kms 1 for the OB supergiants, <3kms 1 for the two cool hypergiants and <10kms 1 for the 
Wolf-Rayet W239. 

'The broad, weak Paschen-series absorption lines in W30a precluded accurate RV measurement, while RV measurements for the 
eclipsing binary W13 discussed in Section [3. ll will appear in a forthcoming paper. 
c The line(s) used for RV measurement. 

d The range of measured radial velocities is given for targets in which statistically-significant variations are detected. 

e Displays a distorted Pa-16 line due to blending with C IIIA8500 and a weak Pa-15 line that yield discrepant radial velocities 

compared to the other Paschen-series lines. 

^ The C III A8500 and C III A8664 lines seen in the spectrum of the Wolf-Rayet W239 give discrepant RV measurements 
(~25kms _1 ; see discussion in Section f3.2p . The value for C III A8664 is given here. 

9 The cool hypergiants W243 and W265 display rich absorption-line spectra, with a large number of potential lines suitable for RV 
measurement. The values given here are the average of nine of the strongest unblended N I absorption lines from multiplets 1 and 8 
(Moore 1945). Note that in both objects species-to-species variation in RV is seen, with clear evidence of velocity stratification 
in the photosphere of W265 (see discussion in text). 

The range of observed radial velocities in W265 varies depending on the absorption line used to make the measurement. See 
Section 



